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Abstract 

In 1990, Kolesova, Lam and Thiel determined the 283,657 main classes 
of Latin squares of order 8. Using techniques to determine relevant Latin 
trades and integer programming, we examine representatives of each of 
these main classes and determine that none can contain a uniquely com- 
pletable set of size less than 16. In three of these main classes, the use of 
trades which contain less than or equal to three rows, columns, or elements 
does not suffice to determine this fact. We closely examine properties of 
representatives of these three main classes. Writing the main result in 
Nelder's notation for critical sets, we prove that scs(8)=16. 



1 Introduction 

A Latin square L of order n is an n x n array of entries {(i, j; k)} such that each 
row and column of L contains each of n possible elements exactly once. We 
will refer to the Latin square based on the group table of Z„ simply as 1 n . A 
uniquely completable (UC) set U is a subset of a Latin square L such that L is 
the only superset of C which is a Latin square. A critical set C of L is a subset 
of L such that C is uniquely completable and no subset of C has this property. 
Critical sets were introduced by Nelder |17|. 

A related concept is that of the Latin trade. In this paper, we will consider 
a Latin trade to be the set of entries in which two Latin squares of the same 
order differ. An intercalate is a Latin trade of size 4. The connection between 
uniquely completable sets and Latin trades is well-known, for example see |12j . 
8 . It is expressed in the following lemma. 

Lemma 1 In any uniquely completable set U for a Latin square L, each trade 
in L must intersect U in at least one position. 
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scs(n) is a function denned by Nelder as the smallest size of a critical set 
in any n x n Latin square. Nelder determined that scs(3)=2, Curran and van 
Rees [B] that scs(4)=4 and scs(5)=6, Bate and van Rees [2] that scs(6)=9, and 
Adams and Khodkar P that scs(7)=12. Nelder [TQ, Bate and van Rees [2], and 
Mahmoodian JS] independently conjectured that scs(n) = \Jj-\- Here we show 
that scs(8)=16. 

For a set S in a Latin square L we define sets for each row i, column j and 
element k. Let Ri(S) = {k\ (i,j;k) G S}, C 3 {S) = {k\ (i,j;k) 6 S}, and 
E k (S) = {(i,j)\ (i,j;k) E S}. So R t (S) (Cj-(S)) is the set of elements which 
appear in row i (column j) of S and Ek(S) is the set of positions where the 
element fc appears in S. 

Now define R{S) = \{i\Ri( s ) + c ( s ) = \{j\Cj(S) # 0}|, E(S) = 
\{k\E k (S)?9}\. 

We will write that a trade T has R(T) rows, C(T) columns and E(T) ele- 
ments. 

Based on results by Cavenagh 0, and considering the paper by Horak, 
Fleischner, and Aldred JUj, the author decided to conjecture the following. 

Conjecture 1 In a subset S of L, a Latin square of order n > 3, S < \J^\, 
there exists a a Latin trade T 6 L such that T OS = and R{T) < 3, C(T) < 3, 
or E{T) < 3. 

2 Integer programming methods 

A survey of defining sets by Donovan et al (E| includes a description of an algo- 
rithm called "Algorithm B" and the minimising technique reproduced verbatim 
here. Smallest defining sets for designs are a concept analogous to smallest crit- 
ical sets for Latin squares, and blocks are here analogous to entries in Latin 
squares. 

Algorithm 1 (B) Find some trades in D and put these in a list T . Now 
perform steps (1) and (2) repeatedly until S has only one completion. 

(1) Form the integer programme corresponding to T , find an optimal solu- 
tion, and form the set of blocks S corresponding to this solution. 

(2) If S does not complete uniquely, then completions not equal to D define 
additional trades in D, and these are added to T. 

The running time for step (1) can be improved by minimising T. That is, 
if T a 6 T and Tf, £ T , and if T a C Tf,, then Tb can be removed from T . 

This algorithm is suitable for finding a smallest defining set for a design, but 
it is not suitable for finding the size of the smallest critical set in Latin squares 
of order 8, as it requires too much computation time. In this paper, we are 
concerned only with the size of the smallest critical set. 

In the paper by Donovan et al it is stated that "Algorithm B yields only a 
single smallest defining set". However, it can be easily extended with a third 
step, as follows, to find all smallest defining sets. 
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(3) If S has unique completion and contains s blocks, add an additional 
constraint C to the integer programme which ensures that less than s of the 
blocks in S can occur. Go back to step (1). 

Similarly to Algorithm B, if T is the complete set of Latin trades in a Latin 
square L, then the optimal solution to the following integer programme is the 
size of the minimal critical set in L. 

Minimize: ^2 xeL C x 
Subject to: 

for each reT^^C.M 

where C x is 1 if x € C and otherwise. 

The application of Conjecture 1 means that in order to determine the value 
of scs(n) for small values of n, we should first find representatives of each main 
class (JUL DP) 01 Latin square of order n and for each square, find the Latin 
trades on 3 rows, columns and elements and write a 0-1 integer programme 
as above. However, for our purposes it is not necessary to calculate the exact 
value of the smallest critical set for every main class. In general, it was found 
that limiting the size of the trades to a value determined by trial and error and 
considering trades only with less than or equal to 3 rows, columns, or elements 
was effective. Thus the following algorithm was used for a Latin square L of 
order n. 

Algorithm 2 Find all trades with less than or equal to 3 rows, columns, or 
elements. This is achieved by considering all completions of L \ S, where S is 
a set of 3 rows, columns, or elements of L. When S ranges over all possible 
sets of 3 rows, columns, or elements of L, then we can find all such trades by 
considering the difference between L and each possible completion. See Bean ^ 
for completion methods. Minimize the list of trades as in Algorithm B. 

In the following, the initial value of z and the increment are chosen to min- 
imize running time. 

(1) Create an integer programme IP as above using the trades of size less 
than z. (Adding an extra constraint to ensure that the sum of the C' x s does not 
equal or exceed \_^-\ speeds up the process.) Solve IP. 

2 

(2) If IP has a solution with size < [^J then increase z and repeat step 1, 
until z = 3n. Otherwise stop. 

(3) Similarly to the new step (3) above, check the solution S of size s < [ 7 -^-\ 
for unique completion. If it does not have unique completion, remove it from 
future consideration by adding a constraint which ensures that less than s of the 
cells of S are considered in future iterations. 

We used CPLEX or BonsaiG |2J to solve the integer programme. 

scs(n) may be determined by considering all the main classes of Latin squares 
of order n, as Donovan et al [7] showed. Kolesova, Lam and Thiel p|] determined 
the 283,657 main classes of Latin squares of order 8. These, and the main classes 
of smaller orders, were obtained from McKay's website 14 . 
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For n — 6, using Algorithm 2 and limiting the size of the trades to 10 or 
less, we can determine that scs(6)=9 in approximately 5 seconds on an Athlon 
1200Mhz computer. Twelve main classes must be examined, thus each main 
class takes about half a second. For n — 7, by limiting the size of the trades 
to 11 or less, we can determine that scs(7)=12 in 34.5 minutes on the same 
computer. We examine 147 main classes, thus each main class takes about 12 
seconds. 

As might be expected, there is a strong correlation between the number of 
intercalates in a Latin square of order n and the time taken to prove the non- 
existence in the square of a critical set of size less than |_xJ" ^or example, 
there are two main classes of 7 x 7 Latin squares with no intercalates; in the 
computation mentioned above, one of these main classes, Z7, accounts for 13.8 
minutes of the total 34.5 minute computation. The other main class contains 12 
Latin trades of order 6 (the size of the smallest trade in Z7 is 9) and accounts 
for 2.2 minutes. The one main class with 1 intercalate takes 1.1 minutes, and 
every other main class takes 45 seconds or less. In the 6x6 Latin squares, the 
effect is less apparent, but the four main classes with less than [^-\ intercalates 
are the four Latin squares which take the most time. 

For a given order 8 Latin square, it is much easier to test whether a critical 
set of size less than 16 exists than it is to find the size of the smallest critical set 
in that square. This is strictly a constraint satisfaction problem rather than an 
integer programming problem, because the minimization step is not necessary. 
Increasing the number of constraints, for this problem, seems to in general slow 
solution down. We began by limiting the size of the trades to 10, which was 
sufficient for more than 95% of Latin squares of order 8, before increasing the 
limit to 12, 14 and then 24. 

For the first 138,800 main classes (arranged in descending order of the num- 
ber of intercalates) the computers used were an Athlon 1200Mhz, and computers 
with 5 Athlon MP 1667Mhz CPUs at IPM, Tehran. CPLEX was used. After 
this a network at the University of Queensland consisting of 128 Sun Fireblades 
and 33 dual Pentium III 800 MHz computers were used with BonsaiG to com- 
plete the calculation, with completion of some gaps by a Pentium IV 2.4 Ghz 
at IPM. The programmes written can be obtained by emailing the author. The 
first part of computation took about three months, and the second part about 
one month. 

As Conjecture 1 was true for n — 7 it was conjectured that it would be true 
for all n. However, while determining the main result of this paper, the author 
discovered the following three main classes X, Y, and W (see Table 1). In A 
there exist 4 sets X 1 ,X 2 ,X 3 , and X4, each of size 15, such that for each Latin 
trade T in X u 1 < i < 4 such that T n X t = 0, R(T),C(T) and E(T) are all 
greater than 3. In other words, for each of these sets, every trade in the square 
with less than or equal to three rows, columns, or elements intersects this set. 
Similarly, in Y there exist 12 sets of size 15 Y\, . . . , Y12 with this property, and 
in W the set of size 15, W\. These are the only three 8x8 main classes with 
this property, thus step (3) of Algorithm 2 is used only for these main classes. 
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These sets are presented in Tables 2, 3 and 4. 

Thus even if a method could be found to generalize the trades constructed 
by Cavenagh to all trades on three rows, columns, or elements, such a method 
would not be enough to prove that scs(n) = \_\\ ■ 

3 Properties of the exceptional main classes 
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Table 1 : The three main classes X, Y, and W 
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Table 2: Xi,X 2 ,X 3 and X 4 
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Table 3: Yx, . . . ,Y 12 
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Table 4: Wi 



We now examine properties of the main class representatives X, Y and W. 

X and Y are in the class of Latin squares Meynert ^Hj refers to as res — 
symmetric. Thus, k) £ L implies entry (j, k;i) € L and entry (k, 6 L. 
Although they are not presented in this way here, they appear in this class. 

Each of the sets X t and Y % has 11,662,776=2 3 .3 2 .161, 983 completions. W x 
has 7,075,188 = 2 2 .3 5 .29.251 completions. 

Each of the sets X i: Yi, and have R(j),C(j), and E(j) for j = 1, . . . , 8 
equal to 1 or 2. Since the size of each set is 15, exactly one row and column 
contains one entry, and one element occurs once. 

X contains 21 intercalates, Y contains 7 intercalates, and W contains 9 
intercalates. Note that X has three times as many intercalates as Y and Y has 
three times as many sets as X with the property above. 

The sets Yi are symmetric about the main diagonal. 

For both X and Y, there are seven intercalates intersecting the entry (8, 8; 8) 
and in both cases the sets Xi and Yi use in total exactly 43 entries from the 
squares X and Y respectively. In X% each entry present occurs once, except 
(8, 8; 8) which occurs 4 times and a transversal (a set in which every clement, 
row, and column occurs exactly once) in which each entry (except (8, 8; 8)) 
occurs 3 times; in Yi each entry present occurs 4 times, except (8, 8; 8) which 
occurs 12 times. 

X and Y seem to have a similar spectrum of possible critical set sizes - in 
both cases, the spectrum found so far is 21 to 28. This leads naturally to a 
question, given that they share so many other properties. 

Question. Do there exist two main classes of Latin squares of the same 
order such that there is a mapping between critical sets from one main class to 
the other? 

4 Ideas for future research 

Further research could include looking at the 1707 main classes of order 9 which 
contain no intercalates, to check whether any critical sets of size less than [^-J 
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exist there. Alternatively, the computation could be repeated to verify "Con- 
jecture 1" of Bate and van Rees [3] that the critical set of size |_^j-J ex ists only 
in the main class Z n , for n = 8. 
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